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We have investigated the effect of epitaxial strain on the magnetic properties and B-site cation ordering in 
multiferroic Bi2FeCr06 (001) thin films using a density-functional theory approach. We find that in thin films 
with rock-salt ordering of Fe and Cr the ground state is characterised by C-type anti-ferromagnetic (AFM) 
order. This is in contrast to the bulk form of the material which was predicted to be a ferrimagnet with G- 
type AFM order. Furthermore, the cation ordered thin-films undergo a transition with epitaxial strain from C 
to A-type AFM order. Other magnetic orders appear as thermally accessible excited states. We also find that 
B-site cation disordered structures are more stable in coherent epitaxial strains thereby explaining the lowered 
magnetic moments observed in these samples at room temperature. Strain varies both the sign as well as strength 
of the Fe-Cr superexchange coupling resulting in a very interesting phase diagram for Bi2FeCr06 thin films. 

PACS numbers: 77.55.Nv, 75.50.Gg, 75.50.-i, 75.10.-b, 75.25.-j, 75.30.-m, 75.80.+q 


Multiferroic materials are characterised by coexistence of 
ferroic orders such as ferroelectricity, ferromagnetism and/or 
ferroelasticity along with a coupling of at least two of these or¬ 
ders. These properties have enormous technological implica¬ 
tions QUI. Hence, much attention has recently been directed 
towards understanding the origins of multiferroicity and the 
design of multiferroic materials 130. 

Presently there are only a few examples of candidate mul- 
tiferroics most of which exhibit weak magneto-electric cou¬ 
plings, particularly at room temperature. The desire for a mul¬ 
tiferroic with significant values of both polarisation and mag¬ 
netisation at room temperature motivated the ah initio design 
of the double-perovskite Bi2FeCr06 (BFCO) derived from the 
parent multiferroic BiFeOg (BFO) 0191. Both compounds 
posses rhombohedral distorted structure in bulk form with 
space group R?>c. The structure consists of Fe and Cr cations 
occupying the 5-sites on alternate (111) layers (rock-salt or 
double-perovskite ordering). A G-type anti-ferromagnetic 
(AFM) order was predicted in bulk BFCO with Fe and Cr sub¬ 
lattices having opposite spins. This mirrors the ordering seen 
in the parent compound. The source of ferroelectric property 
of both BFO and BFCO is the lone pair on Bi whereas the 
ferromagnetism arises from the superexchange interaction of 
magnetic ion pairs Fe((i^)-Fe((i^) and ¥Q{d^)-Cr{d^), respec¬ 
tively O . In BFCO this implies a bulk magnetic moment of 
2/iB per formula unit making it a ferrimagnet. 

Experimental efforts to realise BFCO have resulted in sev¬ 
eral detailed studies of epitaxially grown BFCO on oxide sub¬ 
strates nMEi- These studies have confirmed the existence 
of a room temperature magnetisation as well as large polari¬ 
sation. However, the magnetic moments measured are gener¬ 
ally up to an order of magnitude lower than the theoretically 
predicted value in most films whereas the polarisation remains 
relatively unaffected. The electric and magnetic properties of 
BFCO films have been studied as function of epitaxial strain 
and orientation. It has emerged that under most conditions the 
double-perovskite structure predicted by baettig et al. (H is 
only metastable in the films which instead have the transition 


metal ions disordered over the 5-sites. This disordering of 
Fe and Cr ions has been suggested to be responsible for the 
low magnetisation seen due to local cancelation of moments. 
Recently, Khare et al. C3 observed that the magnetisation in 
epitaxially grown Fa-doped BFCO (001) thin films at room 
temperature is not only significantly reduced (< 0.3/iB/f-ti.) 
compared to the bulk value but increases in magnitude from 
tensile to compressive strains. The connection between the 
magnetic moments in the films, epitaxial strain and the 5- 
site cation ordering is as yet unclear and can be best pro¬ 
vided by first-principles simulations. However, no such stud¬ 
ies presently exist in the literature. 

In this study, we investigate the origins of the low mag¬ 
netic moments and the strain dependence of the moments in 
epitaxially grown BFCO (001) thin films by employing first- 
principles density-functional theory (DFT) simulations. By 
considering a variety of magnetic and 5-site cation order¬ 
ings in the thin films we establish that contrary to the bulk 
the ground-state, AFM order in these films is not of the G- 
type. Throughout coherent strains we observe the C-type 
AFM ordering to be the most stable while other orders ap¬ 
pear as thermally accessible excited states. Furthermore, we 
also find that structures with anti-site defects (see Fig[^ are 
more stable than the rock-salt ordered structure through all 
coherent strains. For the rock-salt ordered structure we find 
that epitaxial strain not only induces a transition in the mag¬ 
netic order from C-type AFM to A-type AFM but also lowers 
the gap to the G-type AFM and ferromagnetic (FM) orders. 
Since the latter two are associated with non-zero moments 
this means a potential enhancement in the magnetic moment 
with strain at room temperature. The drastic variation of the 
energies of these magnetic orders can be understood by the 
variation of the in-plane and out-of-plane exchange coupling 
constants. This is possible because the superexchange inter¬ 
action in BFCO involves sl d^ — d^ pair which, depending 
on the Fe-O-Cr bond angle, can be either FM or AFM in na¬ 
ture OMSl- We have confirmed this variation by extracting 
the coupling constants from supercell calculations at differ- 
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FIG. 1. (a) Different types of super cells considered in the calculations differing in : (i) magnetic ordering, with up and down arrows 
representing up and down spins, respectively, and (ii) cation ordering giving DO (double perovskite), D1 and D2 structures., (b) Energy vs 
strain for all magnetic and chemically ordered structures of BFCO. Energies are plotted relative to the DO C-type structures (broken line) for 
the DO structure (bottom panel), D1 structure (middle panel) and the D2 structure (top panel). In all panels the following symbols are used for 
the magnetic orders : asterisks for C-type, diamonds for G-type, circles for A-type and squares are for ferromagnetic order. 


ent strains. We have also used the coupling constants thus 
extracted in a finite temperature 3D Monte Carlo simulation 
to calculate the layer distribution of the thin films over vari¬ 
ous magnetic orders. Thus we can estimate the magnetisation 
of the sample at room temperature under coherent epitaxial 
strain averaged over different cation orderings (Fig[^). We 
find the estimated magnetisation to be of similar magnitudes 
as measured in experiments HU. However, our results indi¬ 
cate almost no strain-dependence of the moments at least in 
the coherent strain region. Any dependence must either result 
from the differing proportions of cation ordering in the ex¬ 
perimental samples or from spin canting which has not been 
accounted for in this work. 

Experimentally, it was found that BFCO (001) films grow 
with a tetragonal (Pbnm) structure ifTSl IVfl under coherent 
strains. Thus, for our calculations we constructed a 20-atom 
^/2 X ^/2 X 2 tetragonal supercell (inset in Fig[^), starting 
from a simple cubic double perovskite structure, to allow for 
appropriate magnetic ordering of ions along (111) direction. 
Using the experimental in-plane pseudo-cubic lattice param¬ 
eter a cub = 3.93 A ITSJIIII as reference we generated struc¬ 
tures mimicking the epitaxially-strained films by varying the 
in-plane lattice parameter. The in-plane lattice parameters for 
our supercell were set to a = 6 = ^/2 x acub while the c 
parameter was relaxed for each in-plane strain. The epitaxial 
strain was defined as e = (a — aref)lb.ref where aref is the 
unstrained lattice parameter. 

Our calculations employed a spin-polarized GGA-\-U 03 - 
IID approach using the Perdew-Burke-Ernzerhof 1^ 1^ 
exchange-correlation functional within the framework of the 
Quantum-ESPRESSO code 1^ . We chose U = 3.0 eV 
which reproduced well the band structure obtained in previous 


LSDA-fU calculations on BFCO (111251. For further details 
we refer the reader to the supplementary information (SI). 


We first investigated the effect of epitaxial strain on the G- 
type AFM ordered BFCO since this was the order predicted 
for the bulk. For each value of the strain e the atomic po¬ 
sitions as well as the out-of-plane lattice parameter were re¬ 
laxed yielding a value c^^^(e) that minimises the total energy. 
The calculated (Fig 3(a) in SI) are within 1.0% of the 

experimentally reported values for similar strains. Simi¬ 
lar agreements have been obtained in previous first-principles 
studies on epitaxially strained films (iSh giving us confidence 
in the method employed presently. We note here that in the 
experiments ca coherently strained films formed only in a 
window of 4% about the reference strain. Hence our results 
would be most relevant in this strain window. Our calculations 
also extended beyond this window in order to clearly illustrate 
the effect strain has on the various interactions in the system. 


The average magnetic moment per formula unit (f.u.) was 
calculated at each Cmin- But the spontaneous magnetisation 
remained fiat throughout the coherent strain region at the bulk 
value of 2fiB- This trend clearly disagrees with the observa¬ 
tion of lowered moments in experiments CD and motivated us 
to consider the possibility of different ground state magnetic 
and cation ordering. 

We considered three possible types of AFM ordered struc¬ 
tures - A, C and G-types - and a FM ordered structure (see 
Figjjb)). In each of these magnetic orders we further allowed 
for cation disorder by considering two different arrangements 
of the Fe/Cr ions in the supercell (anti-site defects) besides the 
double-perovskite structure, DO. These two different types of 
defect structures (Figj^b)) are referred to as D1 and D2 struc¬ 
tures below. These defects differ in the occupation of the per- 
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ovskite 5-sites in the supercell. Note that only a few of these 
structures result in a non-zero magnetic moment, viz. G-type 
in DO, A-type in D1 and C-type in D2. Needless to say, in 
all structures the FM arrangement leads to a non-zero mag¬ 
netic moment as well (8 /i^/f.u.). Experimentally grown thin 
films may present lower concentration of anti-site defects than 
those considered here, implying the need for larger super-cells 
in calculations. However, in the spirit of keeping the analysis 
simple we restrict our studies to a very limited set of disor¬ 
dered structures within the 20-atom super-cell. 

Taking the DO (double-perovskite) supercell we considered 
each one of the magnetic orders in Figjjb) and carried out 
relaxations similar to the one described above. The min¬ 
imised energy at each epitaxial strain is plotted against strain 
in Fig 1 1(b) I Interestingly, the C-type AFM order emerged as 
the ground state from -4% strain upto 1.5% strain. This is 
in contrast to what is anticipated from bulk BFCO. The G- 
type and A-type AFM, and the FM orders appear as excita¬ 
tions from this phase. In particular, the G-type is very close 
in energy to the ground-state order being easily thermally ac¬ 
cessible. From 1.5% to 4% strain the A-type AFM order has 
the lowest in energy with the G-type AFM and FM ordered 
excitations being nearly degenerate. This clearly indicates the 
possibility of a strain-induced transition in the magnetic or¬ 
der in BFCO thin films. Compressive strains seem to pre¬ 
fer the C-type AFM order while tensile strains induce A-type 
AFM order. These results are interpreted in terms of the super¬ 
exchange coupling constants below. 

Similar calculations were also carried out for D1 and D2 
structures. The energy versus strain plots of all structures 
(Fig |l(b)| ) suggest an interesting phase diagram of different 
possible orders. It was found that under all strains the D1 
structure in the C-type AFM ordering was the most stable. 
The ground-state order in the D1 structure yielded a zero mag¬ 
netic moment and dominates the phase diagram in the exper¬ 
imental range of strains. This could explain the lowered ob¬ 
served moments in the films compared to the bulk value. In¬ 
terestingly, the ground states for all BFCO structures have a 
zero moment indicating that the experimentally observed mo¬ 
ments arise primarily from thermal excitations. For the anal¬ 
ysis below we focus only on the regions of strains |e| < 2 as 
coherently strained films form only in this window na. 

To explain the strain-dependence of the magnetic moment 
we mapped our 20-atom supercell to a spin-lattice model with 
only the magnetic ions. The unit cell has 2 Fe/Cr layers, la¬ 
beled I and II, and each layer having two (5-) sites, labeled 
a and p. The following Heisenberg hamiltonian describes the 
spin-spin interactions in the lattice 


I qll,a r I , qll,^ , qll,^ 


+ ^i,tk • + ^i!j’,k+i}] 


'II,0 


( 1 ) 


where Jy and J± are the in-plane and out-of-plane coupling 
constants, respectively; i^j^k are site indices and S are spin 
vectors. Given the large magnitudes of the spin we are deal¬ 
ing with = 5/2 and = 3/2) we treat the classical 
version of this hamiltonian choosing the spin vectors to be 
parallel to the z-axis. Using Eq.j^for different magnetically 
ordered supercells it is straightforward extract the values of 
J|| and J± in the DO, D1 and D2 structures in terms of the 
corresponding supercell total energies. The variation of these 


coupling constants with epitaxial strain is plotted in Fig |2(a) 
for the three structures. 


In DO (bottom panel) Jy is positive in the strain window 
—3.8% to 1%, indicating an in-plane AFM coupling between 
Fe and Cr. Simultaneously, in the same strain window, J± 
is negative resulting in an FM coupling between the out-of- 
plane Fe-Cr pair. Thus, the ground state becomes C-type 
AFM in DO for compressive strains. Beyond 1% the coupling 
constants change signs, i.e the J_\_ become positive giving an 
AFM coupling and the Jy becomes negative which in turn 
leads to an in-plane FM coupling between Fe-Cr pair. As a 
result of which the DO system undergoes a C-type to A-type 
AFM transition as seen in the phase diagram above. 

The middle panel of Fig |2(^ shows the variation of the cou¬ 
pling constants in the D1 structure. Here, Jy remains positive 
over all the strain maintaining an AFM coupling between in¬ 
plane Fe-Fe pair as required by the Goodenough-Kanamori 
rules (GK) (271, while at the same time Jj_ remains com¬ 
pletely dominating in the negative region over all kind of 
strain. Hence, an FM coupling between the out-of-plane Fe- 
Cr pair comes into play. Due to the above coupling the ground 
state structure is a purely C-type under compressive as well as 
tensile strain region. 

The top panel of Fig |2(a^ is for D2 where the in-plane order 
is AFM from -3.8% upto -1% strain (Jy > 0) and FM be¬ 
yond that strain ( Jy < 0). However, Jj_ remains positive over 
the whole strain region yielding an out-of-plane AFM order. 
Therefore, there is a G-type to A-type magnetic ground state 
transition in the phase diagram (Fig |l(b)| ) for D2. Since the 
magnititude of Jy of D1 remains larger than DO and D2 the 
ground state is a C-type D1-ordered structure. 

In general, the superexchange interaction between like ions 
(Fe-Fe and Cr-Cr) was seen to be AFM under all strains, as 
dictated by GK. However, the results above indicate that in 
the case of Fe-Cr superexchange pairs there is a strong depen¬ 
dence on the strain, and more directly, on the Fe-O-Cr bond 
angle. It was observed that (Fig|^ in SI) for Fe-O-Cr bond 
angles above ~ 149° the in-plane coupling is FM while for 
lower angles an AFM coupling is preferred. The out-of-plane 
coupling also changes sign but at ^ 147°. Similar sign rever¬ 
sal of the Fe-Cr coupling with increasing bond angles has been 
reported for Fe-doped rare-earth orthochromites earlier from 
Mossbauer effect measurements iMl. An ambiguity in apply¬ 
ing GK to pairs arises because of the differing angular 

dependences of the dominant kinetic and potential exchange 
terms in the superexchange. In this case these dominant path- 
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FIG. 2. Strain dependence of (a) the in-plane (IP, Jy) and out-of- 
plane (OP, Ji_) coupling constants in the DO, D1 and D2 structures, 
and (b) the fraction of magnetic domains obtained from Monte Carlo 
simulations and the magnetisation averaged over all three structures, 
M. (Labels are indicated in the respective figures.) 


ways are Fe(t2^)-0(p7^)-Cr(t2^) for kinetic (AFM) and Fe(e^)- 
0(Pa)-CT(t2g) for potential (FM) exchange. The latter is gen¬ 
erally weaker in strength than the AFM pathway ca. How¬ 
ever, given the difference in strength and orientation of the 
corresponding orbital overlaps (a vv. tt) in the two cases there 
is a strong dependence of the dominant pathway on the ge¬ 
ometry of the bond. Since compressive strain forces lower 
bond angles (^ 149°) for in-plane pairs and larger ones for 
out-of-plane pairs generally the AFM dominates in-plane and 
FM dominates out-of-plane (justifying the stability of C-type 
AFM). Since the FM interactions are rather weak the lowest 
excitation under compressive strains often turns out to be out- 
of-plane AFM (G-type or A-type). For tensile strains the sit¬ 
uation is reversed yielding a preference towards A-type AFM 
ground-state in both DO and D2 with in-plane Fe-Cr interac¬ 
tions. 

The phase diagram in Fig | 1(b) [ indicates that magnetic ex¬ 
citations are thermally accessible in all structures. This mo¬ 
tivated us to study the finite temperature behaviour of the 


thin films. Using Monte Carlo (MC) simulations on a 60- 
layer thick film we sampled the configurations allowed by the 
hamiltonian above at 300K at each strain. For each type of ar¬ 
rangement of the cations we monitored the population of the 
domains yielding non-zero magnetic moment. This is plotted 
as the fracti on of the layers belonging to a domain versus the 
strain in Fig 2(b) along with the averaged magnetisation (M). 
The drastic reduction in the magnitude of the magnetisation in 
the coherent strain window observed in experiments is repro¬ 
duced. Note that the estimated magnitudes are an upper bound 
to the magnetisation assuming alignment among the magnetic 
domains. Nevertheless, the correlation between the trend seen 
in Fig 2(b) and the phase diagram in Fig |l(b)| strongly indi¬ 
cate that our finite temperature predictions consistent with our 
first-principles calculations. 

Our results caution against the direct extrapolation of bulk 
results to epitaxially grown thin films. Thin-film BFCO 
may not have the G-type AFM ground-state expected in its 
bulk form. In contrast, under all coherent epitaxial strains 
BFCO favours a C-type AFM magnetic ordering. Further¬ 
more, the double-perovskite structure is unstable under all co¬ 
herent strains towards cation disordering. In all cases con¬ 
sidered the ground-state has zero magnetisation and magnetic 
moment arises from thermal excitations to the G-type AFM 
and FM orders. We confirm this trend using MC simula¬ 
tions which also yield a value for the magnetisation which 
is significantly lower that the expected bulk value. Our pre¬ 
dictions match quite well with similar observations made in 
recent experiments (UlIISl. In fact, a recent work 1^ mea¬ 
suring XMCD and chemically sensitive hysteresis loops on 
BFCO thin film samples have indeed confirmed the presence 
of both very low magnetic moments as well as cation disor¬ 
der. By analysing the dependence of the exchange coupling 
constants on the strain and bond geometries we conclude that 
the strength and sign of the Fe-Cr interaction is tunable by 
strain suggesting interesting possibilities for BFCO and other 
materials of this kind. 

PLD growth conditions could affect the proportions to 
which the film incorporates the various cation orderings dur¬ 
ing formation. This would result in differences in observed 
magnetisation dependences. However, our work suggests that 
intrinsically BFCO displays a very interesting evolution of 
magnetic and chemical order with epitaxial strain. We hope 
that more consequences to this interesting phase diagram will 
be experimentally probed to reveal other exciting phenomena 
in the material. 
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cussions. The authors gratefully acknowledge USER Bhopal 
for computational resources and funding. PCR would like to 
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Supplementary Information 

DETAILS OF DFT CALCULATIONS 


The spin-polarised DFT calculations in this work employed 
GGA-\-U approach using the Perdew-Burke-Ernzerhof 12^ 
[23l exchange-correlation functional. The simplified version 
of DFT-^f/ functional given by Coccocioni et al (2011211 and 
as implemented in the Quantum-ESPRESSO code (^ was 
used. We chose a value of 3.0 eV for the Hubbard parameter U 
which reproduced well the band structure obtained in previous 
LSDA-^f/ calculations on BFCO (811251. But, the Hubbard U 
parameter can, in principle, be calculated ah initio using the 
linear response formalism of Ref. (20l 1^ . We found the U 
values to be of 4.25 eV (Fe) and 3.07 eV (Cr) in BFCO. How¬ 
ever, these neither altered the band structure significantly nor 
affected the phase diagram for the strained thin films beyond 
the accuracy of our calculations. So, to stay consistent with 
literature, we proceeded with U = 3 eV. We used a 8 x 8 x 8 
Monkhorst-Pack k-point mesh for Brillouin zone integration 
and a plane-wave basis with a 58 Ry kinetic-energy cutoff 
along with a 600 Ry cutoff for the charge-density. Ionic cores 
were modelled by ultra-soft pseudopotentials. Structures were 
relaxed until the forces are less than 13 meV/A. 


STRAIN-DEPENDENCE OF OUT-OF-PLANE LATTICE 
PARAMETER 


The calculated Cmin (e) versus in-plane lattice parameter (a) 
along with the experimental result Khare et al ca is plotted 
in Fig. 3(a) within the coherently strained region. At the com¬ 
pressive strain the calculated out-of-palne lattice parameter 
(Cmin(e)) is 1.01% larger than the experimental value and at 
tensile strain it underestimates 0.56% to its corresponding ex¬ 
perimental value but follows the experimental observed trend. 
Similar errors have also been reported in previous GGA stud¬ 
ies on similar materials. See, for instance, the work by Dobin 
et al (26l . The inset in Fig. |3(a)| is the structure of the BFCO 
thin film unit cell where the B site cations are maintained al¬ 


ternatively to get the required magnetic ordering. Fig. |3(b) 
shows the strain dependence of Cjnin foi* all magnetic types 
considered for the double perovskite (DO) structure. There is 
very little effect of the change in magnetic order on the lattice 
parameters. 
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(b) 

FIG. 3. Variation of out-of-plane lattice parameter with the in-plane 
lattice parameter. The curve with circlular symbols represents the 
GGA+U calculated out-of-plane and the line with the square points 
denotes the experimental one.The inset double perovskite structure 
is the unit cell used for the calculations 


STRAIN-DEPENDENCE OF THE Fe - O - Cr BOND ANGLE 


The Fe-O-Cr bond angle was monitored at every strain for 
all magnetic orders. Fig.Qa) shows this dependence for the 
DO structure. As the ground-state changes from C-type to A- 
type AFM (around 1% strain) the in-plane angle increases to 
^ 150° while the out-of-plane angle is lowered to ^ 146°. 


Strain dependence of bond-angles are also shown for D1 
and D2. In D1 the in-plane and in D2 the out-of-plane an¬ 
gles are for Fe-O-Fe or Cr-O-Cr, i.e. between like cations. 
These do not undergo any change in sign with strain and re¬ 
main AFM consistent with the Goodneough-Kanamori rules. 
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(b) 


(b) 


FIG. 4. (a) Variation of out-of-plane (OP) and in-plane (IP) Fe-O-Cr 
angle under epitaxial strain for different types of magnetic orders in 
DO structure. The vertical line marks the transition strain where the 
film goes from C-type to A-type AFM order, and (b) Variation of OP 
an IP coupling constants with the corresponding Fe-O-Cr bond angle. 
Note that the direction of increasing strain is along the direction of 
increasing IP angle and decreasing OP angle. 


MONTE CARLO SIMULATIONS 


The following Heisenberg hamiltonian describes the spin- 
spin interactions in the lattice 


I qll,a r I , qll,^ , qll,^ 

^ ^ ^i-l,j,k ^i,j-l,k ^ ^i-l,j-l,kU 

+ E + ^ij,k+l} 


+ ^i,tk • + ^i!j’,k+i}] 


'II,0 


( 2 ) 


FIG. 5. Variation of out-of-plane (Fe-O-Cr) and in-plane (Fe-O-Fe or 
Cr-O-Cr) angle under epitaxial strain for different types of magnetic 
orders in (a) D1 structure, and (b) D2 structure. 


where Jy and J± are the in-plane and out-of-plane coupling 
constants, respectively; i^j^k are site indices and S are spin 
vectors. Given the large magnitudes of the spin we are deal¬ 
ing with (cr^® = bl2 and = 3/2) we treat the classical 
version of this hamiltonian choosing the spin vectors to be 
parallel to the 2 ;-axis. Using Eq.j^for different magnetically 
ordered supercells it is straightforward extract the values of 
J|| and Jx in the DO, D1 and D2 structures in terms of the 
corresponding supercell total energies. For instance, for the 
DO structure we have 


Ea — Eg Eg — Eg 


(3) 


Similar expressions can be obtained for the coupling constants 
in the D1 and D2 structures as well. 

Using Monte Carlo (MC) simulations on a 60-layer thick 
film we sampled the configurations allowed by the hamil¬ 
tonian above at 300K at each strain. Furthermore, the MC 
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moves were restricted to those that flip spins of an entire sub¬ 
lattice in a given layer. This was done in order to accelerate the 
sampling since we were primarily interested in the population 
of domains of each magnetic order in the thin-fllm. For each 
type of arrangement of the cations we monitored the popula¬ 
tion of the domains yielding non-zero magnetic moment. The 
fraction of the layers belonging to a domain versus the strain 
is plotted in the main text in Fig.2(b). 

The average magnetic moment can be estimated by the sum 
of the contributions of the magnetic domains from each struc¬ 
ture weighted by the probability of the thin-fllm occurring in 
a given structure 

M(e) = J2w‘'{e) (4) 

where n\{e) is the fraction of X-type magnetic domain in u- 
type structure at the strain e and is the magnetic moment 
per formula unit associated with a X-type magnetic domain in 
z/-type structure. X and v run over the different magnetic and 
structure types DO, D1 and D2 (see main text). w^{e) is the 
weight associated with the occurrence of the z/-type structure. 
The latter is approximated as the probability of occurrence of 
the ground-state magnetic domain in that structure 

Y,exp{-/3E’^) 

w’'ie) « ^^-; Z = Y,Y.^xp{-^E^,) (5) 

V T 

where refers to the energy in the magnetic domain r (r = 
A,C,G & Ferro) of the z/-type structure. The denominator in 
the above equation is chosen to ensure normalisation. M is 
also plotted as a function of strain in Fig.2(b) of the main text. 

Note that the estimated magnitudes are at most an upper 
bound to the magnetisation assuming alignment among the 
magnetic domains. The latter is obviously not necessary since 
at flnite temperature out-of-plane fluctuations may lead to de¬ 
struction of any order in the thin-fllm. Also in our choice of 
MC moves we have ignored in-plane fluctuations of the spins 
which might destroy even in-plane order and reduce magneti¬ 
sation further. 
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